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Abstract 

In this paper, we investigate the integrability aspects of a physically important nonlinear oscilla- 
tor which lacks sufficient number of Lie point symmetries but can be integrated by quadrature. We 
explore the hidden symmetry, construct a second integral and derive the general solution of this 
oscillator by employing the recently introduced A-symmetry approach and thereby establish the 
complete integrability of this nonlinear oscillator equation from a group theoretical perspective. 
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I. INTRODUCTION 



In this paper, we consider the nonlinear oscillator equation, 

A(t, x, x, x) = (1 + kx 2 )x — kxx 2 + a 2 x = 0, (1) 

where k and a are arbitrary parameters and overdot denotes differentiation with respect 
to t, which was introduced by Mathews and Lakshmanani. Eq.([T]) is obtainable from the 
Lagrangian L = ~ ( '^i^jt ^ which can be considered as the one dimensional analogue of the 

Lagrangian density L = ~ ( 1+ fc^,2 j (d^d^cf) — a 2 <p 2 ) appearing as a nonpolynomial model in 
quantum field theory^. Here a and k are parameters. Eq.fll]) admits the general solution of 
the form x(t) = A sin(Qt + <fi) , where A and 4> are two arbitrary constants and the frequency 
f2 is related to the amplitude A through the relation Q 2 = ■ While Mathews and 
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Lakshmanan have also proved the quantum exact solvability of Eq.(JI]) , its three dimensional 
version was shown to be solvable both classically and quantum mechanically by Lakshmanan 
and Eswaran^. Later Higgs and Leemon studied the classical and quantum dynamics of the 
oscillator (|T|) on the spherical configuration spaced 

In recent times multifaceted investigations have been carried out on this equation and its 
higher dimensional generalizations at the classical as well as quantum levels. For example, 
Carinena and his collaborators have generalized the above nonlinear oscillator model to 2- 
and 3-dimensions and formulated an n- dimensional integrable analogue of it and unearthed 
several mathematical properties associated with the underlying equations, including general 
solution, Lagrangian and Hamiltonian structures, conserved quantities and so on 6 . The 
authors have also demonstrated that the underlying system is a superintegrable one. A 
geometric interpretation of the higher dimensional system was also proposed in relation 
with the dynamics on spaces of constant curvature 6 -. Interestingly, Eq.([T]) is linearizable to 
harmonic oscillator through a nonlocal transformation. The two dimensional generalization 
is also linearizable^. By making use of the geometric properties of the kinetic energy of the 
corresponding Hamiltonian, the associated quantum version of system (TjQ) was studied and 
considered as a deformation of the linear harmonic oscillator-^. The quantum version of the 
above nonlinear harmonic oscillator has also been shown to be factorizable and its shape 
invar iance property was brought out. 

Very recently Midya and Roy have generalized this nonlinear oscillator and analyzed the 



underlying Schrodinger equation under various aspects including exact solvability, quasi- 
exactly solvability and non- Hermit ian invariants^. The coherent states for this nonlinear 
oscillator has also been constructed and its dynamics was studied in Ref. 
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The studies 



recalled above indicate the necessity to make further progress on the physical and mathe- 
matical structures behind the model ([I]). 

The aim of this paper is to establish the complete integrability of Eq.flT]) through a Lie 
point symmetry analysis which is a nontrivial task since the nonlinear oscillator equation 
admits only the time translational symmetry. One can deduce the first (energy) integral 
from the time translational symmetry which in turn can be integrated to yield the general 
solution. One can also rewrite the first integral as a linear first order equation in some new 
variables, which can be integrated to give the second integral (see Appendix A for details). 
However, we aim to explore the second symmetry which can be used to construct the second 
independent integral and leads to general solution without any quadrature. 

To explore the second symmetry and its associated integral we consider the recently 
introduced A-symmetry method^"—. The method of finding A-symmetries for a second-order 
ordinary differential equation (ODE) has been discussed recently by Muriel and Romero^ 5 -. 
Once the A-symmetries are found out then one can construct the integral directly from the 
A-symmetry by a four step algorithm. The underlying integrating factor can then be deduced 
from the integral just by differentiating the latter with respect to x. It has been shown that 
in certain specific examples 13 that eventhough the determining equations for classical Lie 
symmetries admit only trivial solutions, the corresponding equations for the A-symmetries 
admit nontrivial solutions. On the other hand the function A can be chosen in such a way 
that the corresponding infinitesimals can be calculated. 

As far as the present problem is concerned we find the A functions by solving the A- 
invariant equation. Once the A functions are determined we then proceed to construct the 
integrals through the prescribed algorithm 15 . The resultant integrals provide the general 
solution. 

The plan of the paper is as follows. In Sec. 2, we briefly recall Lie point symmetries 
of Eq.flTJ. In Sec. 3, we present the method of finding A-symmetries and the procedure 
to construct the associated integral from the A-symmetries for a given second-order ODE. 
In Sec. 4, we apply the procedure to Eq.([T]) and obtain two A-symmetries such that the 
associated first integrals are functionally independent. From the latter we deduce the general 
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solution. Finally, we present our conclusions in Sec. 5. 



II. LIE POINT SYMMETRIES OF EQ.® 

In this section we study the Lie point symmetries of Eq.([T]) and show that it admits 
only time translational symmetry. Let the evolution equation be invariant under the one 
parameter Lie group of infinitesimal transformations, 

t = t + e£(t,x) + 0(e 2 ), 

x = x + erj(t,x) + 0(e 2 ), e < 1, (2) 

where £ and rj represent the infinitesimal symmetries associated with the variables t and x 
respectively. The associated infinitesimal generator can be written as 

Eq. ([I]) is invariant under the action of ([3]) iff— 7 ^ 

^ {2) (A)U=o = 0, (4) 

where 

v ® = £+^+* (1) ir+<i w jF (5) 

at ox ox ox 
is the second prolongation in which 

n^ = fj — x£, n^ = r) — x£ — 2x£, (6) 

and dot denotes total differentiation. 

By analyzing Eq.flT]) we get the following determining equations, 

(1 + kx 2 )£ xx + 5kx£ x = 0, (7) 

(1 + kx 2 ){7] xx - 2£ tx ) + 3kx Vx - Akx^t + r)k(l - kx 2 ) = 0, (8) 

(1 + kx 2 ) (2r] tx -£ tt ) + 2kx Vt + 3a 2 x£ x = 0, (9) 

(1 + kx 2 ) 2 r] tt + a 2 x(2£ t - r) x )(l + kx 2 ) - na 2 (l - kx 2 )(l + kx 2 ) 2 = 0, (10) 

where the subscript denotes partial differentiation. Solving Eqs. (!7|)- (|T0|) . we find that the 
only admissible solution is £ = constant and n = 0. As a result one obtains only the time 
translational symmetry, that is 



Consequently one is interested to investigate the existence of more general symmetries to 
establish the integrability of ([!]). 



III. A-SYMMETRIES AND INTEGRALS 

In order to be self-contained, in this section, we briefly recall the notion of A-prolongation 
and the method of finding A-symmetries and constructing the first integral from the A- 
symmetry 1 ^ - — . 

A. A-prolongation 

Let A(t,x,x,x) = be a second-order ODE. A vector field V, is a symmetry of the 
second-order equation if there exists a function such that 

V [x ' (2)] (A(t,x,x,x)) = when A(t,x,x,x) = 0, (12) 

where V' A '' 2 '' is given by 

^ A - (2)] = C(t, x)^ + V iX > m (t, x)-^ + nM 1 * (t, x,x)^ + r,^ (t, x, r.pS. (13) 

with 

V A ' (0)] = v(t,x), (14) 

^,(1)] = (A + X ) v ^(t,x) - (A + \)(£)x, (15) 

^,(2)] = (A + ^,(1)]^^) _ ( A + A )(^ £i (16) 

It has been proved that 14 if V = £(t,x)f + r}(t 

' ~dx ^ ^ 

A- symmetry of the ODE 

x = (J)(t,x,x) for some function A = X(t, x, x) then V — 4- is also a A-symmetry of the 



same equation for the function A = A + —j^p-, where Q = rj — £x is the characteristic of 
V and A = t| + x-^ + (f)(t,x,x)-^r is the vector field associated with the equation under 
consideration 14 . In this case the A determining equation simplifies to 

<f) x + \<J) ± = \ t + x\ x + (j)\ ± + \ 2 . (17) 

In the following we solve the invariant condition (1171) and obtain the A symmetry. 
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We have already observed that Eq.flT]) admits the infinitesimals £ = constant, n = 0. Let 

~ ~ (j) (x, X) 
us consider V = d t is a A-symmetry of (IT]) with A = 0. In this case Q = —x, A = — , 

x 

x( kx^ ct^l 

where x) = — : — r— . Now one can directly check that this A is a solution of Eq. (TT7|) 

(1 + kx 2 ) 

as well. 

If we manage to obtain second solution, A2, which is different from Ai then we can 
construct the second integral from A2 by using the same algorithm given below. In this case 
V = is both a Ai-symmetry and a A2-symmetry^. Since the functions Ai and A2 are 
different the associated integrals are also independent. 



B. First integrals 

The method of generating the first integral from the A- symmetry essentially consists of 
solving the following two equations (see Ref.[15]), that is 

v [Mi)]j = i x + XI A = 0, (18a) 
DJ = 0, (18b) 

where u^'Wl is the first-order A-prolongation of v and D t = ^ + x-^ + (j)(t, x, x) Jr is the total 
differential operator. To begin with let us suppose that w(t, nontrivial integral / 

of t;[ A >W], besides the trivial integral / = t (or a function of t), that is w(t, x, x) is a solution 
of the first-order partial differential equation 

w x + \(t, x, x)w± = 0, (19) 

where subscripts denote partial derivative with respect to that variable. 



Step 1: 

Find a first integral w(t,x, x) which is the solution of the equation (fT9]) . Since there is 
no t derivative in the equation tAWl/ = J x 4- \I ± = 0, the integral of v^- x '^ can be written 
in the most general form as, 

I(t, x, x) = G(t, w{t, x, x)), G w ^ 0, (20) 
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for some function of two variables G(t, w). Now substituting this form in the second condition 
(]18bp . we get 

DJ = G t + D t (w)G w = 0. (21) 

In other words, D t (w) = — q, which will obviously be again a function of the variables 
t and w. We call this function as F(t,w). Consequently, Eq.( l2~Tj) can be expressed as the 
first-order ODE 

G t + F(t,w)G w = 0. (22) 
Note that in the case where w itself is an integral, then D t w = = F(t,w). 

Step 2: 

Evaluate D t w and express D t w in terms of (t, w) as D t w = F(t,w), such that F is only 
a function of t and u. Here D t is the total differential operator. 

Step 3: 

Find a first integral G(t,w) of [d t + F(t,w)d w ]G(t,w) = 0. 
Step 4: 

Then the required integral is I(t, x, x) = G(t, w(t, x, x)). 

Also, fi(t, x, x) = I± is an integrating factor of the given second-order ODE. 

IV. INTEGRABILITY OF EQ. (pQ) 
A. Determination of A-symmetries 

As we mentioned at the end of Sec. 3. 2, A = ^ x : x > is a A-function for ([1]). The explicit 
form of this A is given by 

x ( kx 2 — a 2 



A ' = i W'' (23) 
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To determine a second particular solution of (fl7|) . we assume an ansatz for A in the form 
( a justification for this ansatz is given in Appendix A). 

ai + a 2 x + a 3 x 2 + a 4 x 3 + a 5 x 4 

A — : : ; : — — : — — , (^4) 

Oi + 2 X + 3 X Z + O4X 6 

where cij and fej's, i = 1,2,3,4,5, j = 1,2,3,4, are arbitrary functions of t and x and 
which are to be determined. Substituting this ansatz into Eq.(fl7l) and solving the resultant 
equations we obtain the second independent solution, A2, in the form (the details are given 
in the appendix B). 

. (a 2 - kx 2 )(x + kx 2 tx + kxx 2 + a 2 ktx 3 ) . . 

A2 = • (25) 

(1 + kx 2 ){a 2 x — k(x 3 t + a 2 xtx 2 — a 2 x 3 )) 

Using this second A-symmetry we determine the second integral. We note here that A x 
can also be determined using the above ansatz. Since it has already been found out we do 
not give the explicit form a second time. 

B. Determination of Integrals 

1. Integral from \\ 
Step 1: 

To begin with let us consider the function Ai given in fT23|) . Substituting the latter in 
(19) we get 

x ( kx 2 — CY 2 \ 

w * + i{TTk*r =o - (26) 

We solve the characteristic equation associated with ( 1261) . 

dx dx 

x I kx 2 —a 2 



x \ 1+kx 2 

and obtain a solution w of the form 



(27) 





^• l) = TTfa^' (28) 



s 



Step 2: 

Now evaluating the expression D t w in terms of the variables t and w we find 

D t w = F(t,w) = 0. 



(29) 



Step 3: 

Since F = 0, we find ^=0 which in turn yields G{t, w{x, x)) = G(w(x, x)). 
Step 4: 

Since G turns out to be a function of w only, we simply choose G = w(x,x) and obtain 
the integral as 

= kx 2 -a 2 

1 + kx 2 K ' 

2. Second integral from A2 
Step 1: 

Let us now consider the function A 2 . Substituting (|25|) in (19) we get a complicated 
first-order linear partial differential equation (PDE) of the form 

(a 2 - kx 2 )(x + kx 2 tx + kxx 2 + a 2 ktx 3 ) _ ^ 
x (1 + kx 2 )(a 2 x — k(x 3 t + a 2 xtx 2 — a 2 x 3 )) x 

To construct a solution of the above PDE, we solve the characteristic equation associated 
with (13TI) . namely 

dx (1 + kx 2 )(a 2 x — k(x 3 t + a 2 xtx 2 — a 2 x 3 ))dx 
1 (a 2 — kx 2 )(x + kx 2 tx + kxx 2 + a 2 ktx 3 ) ' 

and obtain a particular solution to the PDE (|3T!) in the form 



a 2 — kx 2 , 



x) = tan" ( ) - yj - - fc ^ 2 t. (33) 
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Step 2: 



Now we evaluate D t (w) in terms of the variables t and w. Here also we find that the total 
differential of the function becomes null, that is D t (w) = 0. This in turn fixes the function 
F(t,w) = 0. 

Step 3: 

As a consequence one essentially obtains a trivial equation in the third step, namely 
G t = 0. 

Step 4: 

Now restricting G(w(x, x))=w(x, x) one obtains an explicit expression for the second 
integral in the form 

a 2 — kx 2 , 



The second integrating factor can be deduced from (IMl) just by differentiating (I34"|) with 
respect to x. 

It is also of interest to note that the second integral I2 given above can be determined 
from the first integral I\ itself, once the later is known explicitly, by the procedure pointed 
out earlier by three of the present authors in Ref.[7]. This procedure also gives a natural 
justification for the ansatz (25). The details are given in Appendix A. 

C. General Solution 

From 7i and I2, we deduce the general solution of (JT|) straightforwardly as 



a 2 



x(t) = Asm(nt + 5), n = \l- — , A=^I U S = I 2 , (35) 

which is the well known solution to ([1]). Note that this integral is related to the standard 
integral, 7i = T^^k , through the relation I\ = k ° ■ 
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V. CONCLUSION 



In this paper we have established the complete integrability of a widely studied nonlinear 
oscillator equation through group theoretical method. The need for this task came from 
the fact that the conventional Lie point symmetry analysis proceduces only one symmetry 
generator. To achieve our goal we have considered the A-symmetry approach. The second 
A-symmetry turns out to be a rational function in x and x. The associated integral also 
turns out to be a complicated rational function. From these two integrals we derived the 
general solution straightforwardly. The integrating factors for the Eq. (TjQ) have also been 
reported. As a by-product of this work we have also shown the utility of the A-symmetry 
approach in solving nonlinear ODEs, where sufficient number of Lie point symmetries do 
not exist. 

Appendix A: A direct method of finding I2 from I\ 

The second integral I2 given in (34) can also be obtained from a knowledge of the integral 
ii given in (30), as follows [18]. 

Let us rewrite the first integral (30) in terms of new variables w and z in the form 

Ikxx at az 

where 

(a 2 -kx 2 ) dt (A2) 

Let us consider a particular solution of (Al) as 

w = hz, (A3) 
where I\ is the integration constant. Now substituting the expression (A2) in (A3), we find 

1 + kx 2 = e hz . (A4) 
Using the first integral (30) we can express x in terms of / and x, that is 

x 2 = (a 2 -h(l + kx 2 ))/k. (A5) 
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Substituting the expressions x and x in the second equation in (A2) and differentiating and 
reexpressing the variables x and x in terms of t and z, we arrive at 

^ = 1((« 2 + I x ) e -^ - h - a 2 e- 2hz ) 1/2 - (A6) 
dt 1\ 

Integrating the above equation, we get 

1 _ x ( Vli(e hz - 1) \ 

Eq.(A7) can be simplified to yield 

'V7i 



J 2 = I 2 ^- h = -Vl.t + tan" 1 {^f-^j 



(A8) 



Substituting the expression (30) for the first integral I\ in (A8) we get the explicit form 
of I2 which coincides exactly with the one given in (34). Once I2 is known the associated 
A-symmetry can be found through the expression [15] 

A 2 = (A9) 

J-2x 

where I2 is given in (34). Substituting the relevant derivatives in (A9) we find that A2 should 
be in rational form. To explore it systematically within the group theoretical framework, 
we assume a rather general form for A2 (vide Eq.(20)) and solve the determining equations 
which exactly leads to the desired expression. 



Appendix B: Method of finding A2 

Substituting the ansatz (!24j) in (ITT)) and rearranging the resultant equation one gets a 
polynomial equation in x with coefficients in Oj and bjS, i = 1,2,3,4,5, j = 1,2,3,4. 
Equating the coefficients of different powers of x to zero one obtains the following set of 
PDEs for a, and b'jS, namely 



(1 + A;x 2 ) 2 (a 5 (a 5 — &4 X )(1 + kx 2 ) + 6 4 (a 5;c (l + kx 2 ) + a 5 kx) — 2b i a 5 kx) 

-b\k(l - k 2 x 4 ) = (Bl) 

(1 + kx 2 ) 2 (a 5 ((2a 4 - b 4t - &3*)(1 + kx 2 ) + 2b^kx) + b 4 (a 5t + a 4x )(l + 
kx 2 ) — 0464^(1 + kx 2 ) + 05^63(1 + kx 2 ) — 2kx{b 4 a 4 + 6305)) 

-26 3 6 4 A;(1 - fc 2 x 4 ) = (B2) 

12 



(1 + kx 2 ) 2 (a 5 ((—b 3t - b 2x + 2a 3 )(l + kx 2 ) - b 4 a 2 x + 3b 2 kx) + a 4 (-b 4t 
—b 3x + «4)(1 + kx 2 ) - a 3 b 4x (l + kx 2 ) + b 4 ((a 4t + a 3x )(l + kx 2 ) - a 3 kx) 
+i> 3 (( a 5t + a 4x)(l + kx 2 ) + a 4 /cx) + b 2 a 5x (l + &^ 2 ) — 2kx[b 4 a 3 + 6304 

+6 2 a 5 )) - (1 - k 2 x A ){b 2 3 k - 6 2 a 2 + 26 2 6 4 A;) = (B3) 

(1 + kx 2 ) 2 (a 5 ((-&2t - &ix + 2a 2 )(l + /cx 2 ) - 2& 3 a 2 x + A^kx) + a 4 ((-& 3i 
-6 2x + 2a 3 )(l + kx 2 ) + 2b 2 kx) + a 3 (-b 4t - 6 3:c )(l + kx 2 ) - a 2 b 4x (l + 
kx 2 ) + b 4 ((a 2x + a 3t )(l + /ex 2 ) - 2a 2 kx) + b 3 (a 4t + a 3x )(l + kx 2 ) + b 2 
(a 4x + a 5t )(l + kx 2 ) + 6ia 52; (l + kx 2 ) - 2kx(b 4 a 2 + b 3 a 3 + b 2 a 4 ) 

+ha 5 ) - 2(1 - k 2 x 4 )(b 2 b 3 k + b x b 4 k - b 3 b 4 a 2 ) = (B4) 

(1 + kx 2 ) 2 (a 5 ((-&it + 2ai)(l + /ex 2 ) - 36 2 a 2 x) + a 4 ((-b 2t - b lx + 2a 2 ) 
(1 + kx 2 ) + 3&iA;x - b 3 a 2 x) + a 3 ((-b 3t - b 2x + a 3 )(l + /cx 2 ) + b 2 kx 
+b 4 a 2 x) + a 2 ((b 3x - b 4t )(l + kx 2 ) - b 3 kx) - axb 4x (l + kx 2 ) + b 4 ((a lx 
+a 2t ){l + kx 2 ) - Zaxkx) + b 3 (a 3t + a 2x )(l + kx 2 ) + b 2 (a 4t + a 3x )(l + 
kx 2 ) + &i(a 5t + a 4x )(l + kx 2 ) - 2kx(b 4 a 1 + b 3 a 2 + b 2 a 3 + 61 a 4 )) 

-(1 - k 2 x%b\k + 2M 3 £; - & 3 « 2 - 2b 2 b 4 a 2 ) = (B5) 

(1 + kx 2 ) 2 (a 4 ((-b lt + 2ai)(l + /cx 2 ) - 26 2 a 2 x) + a 3 ((-6 2t - &ix + 2a 2 ) 
(1 + A;x 2 ) + 2^kx) + a 2 ((-b 3t - b 2x )(l + A;x 2 ) + 2& 4 a 2 x) + ai ((-&&,. 
-b 4t )(l + kx 2 ) - 2b 3 kx) + & 4 a u (l + kx 2 ) + b 3 {a 2t + a lx )(l + /cx 2 ) 
+b 2 (a 3t + a 2x )(l + kx 2 ) + bi((a 4t + a 3x )(l + kx 2 ) — 4a 5 a 2 x) — 2kx 
(ha, + b 2 a 2 + b ia3 )) - 2(1 - A; 2 x 4 )(6 1 6 2 A; - b 2 b 3 a 2 - hb 4 a 2 ) = (B6) 

(1 + kx 2 ) 2 (a 1 ((-b 2x - b 3t + 2a 3 )(l + kx 2 ) + 3b 4 a 2 x - b 2 kx) + a 2 ((-b lx 
—b 2t + a 2 )(l + kx 2 ) + b 3 a 2 x + b x kx) - a 3 b lt (l + kx 2 ) + bi((a 2x + a 3t ) 
(1 + A;x 2 ) - 3a 4 a 2 x) + b 2 ((a lx + a 2t )(l + kx 2 ) - a 3 a 2 x) + b 3 a lt (l + 
kx 2 ) - 2kx(b 2ai + ha 2 )) - (1 - k 2 x A )(kb\ - b\a 2 - 2b x b 3 a 2 ) = (B7) 
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(1 + kx 2 f{a 1 { (-&ix - 6»)(1 + /cx 2 ) + 26 3 a 2 x) + a 2 (2ai - & lt )(l + fcx 2 ) 
+&i((aix + «2t)(l + kx 2 ) - 2a 3 a 2 x) + b 2 a u (l + kx 2 ) - 2kxbi<ii) 

+2(1 - k 2 x 4 )bib 2 a 2 = 



(B8) 



(1 + fcr 2 ) 2 (ai((— bu + ai)(l + A;x 2 ) + &2Ct 2 x) + &i(ait(l + A;x 2 ) — a2« 2 x)) 

+ (1 - Px 4 )6 2 a 2 = (B9) 

To solve these equations we again assume a« and bj's are quartic polynomials in x with co- 
efficients which are unknown functions in t. Substituting these forms in the above equations 
we find each one of the above equations reshapes to a polynomial equation in x. Equating 
the coefficients of different powers of x to zero in each equation one obtains a set of ODEs 
for the differential coefficients. Solving them consistently one finally arrives at the solution 
given in ff25l) . 
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